Mechanisms of enhanced orbital dia- and paramagnetism: 
Apphcation to the Rashba semiconductor BiTel 
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We study the magnetic susceptibility of a layered semiconductor BiTel with giant Rashba spin 
splitting both theoretically and experimentally to explore its orbital magnetism. Apart from the core 
contributions, a large temperature-dependent diamagnetic susceptibility is observed when the Fermi 
energy Bp is near the crossing point of the Rashba spin-split conduction bands at the time-reversal 
symmetry point A. On the other hand, when Ef is below this band crossing the susceptibility 
turns to be paramagnetic. These features are consistent with first-principles calculations, which 
also predict an enhanced orbital magnetic susceptibility with both positive and negative signs as a 
function of Ef due to band (anti)crossings. Based on these observations, we propose two mechanisms 
for the enhanced paramagnetic orbital susceptibility. 

PACS numbers: 75.20.Ck, 71.70.Ej, 85.75.-d, 31.15.A-, 71.15.-m 
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The magnetic susceptibility x is one of the most fun- 
damental physical quantities, revealing invaluable infor- 
mation about the spin and orbital states of solid materi- 
als Q. The sign of x, depending sensitively on the mech- 
anism of the magnetization, plays a crucial role in identi- 
fying the magnetic properties of bulk systems. From this 
viewpoint, solid materials are roughly classified into the 
following three categories ||l[: (i) band insulators with 
weak temperature-independent diamagnetism x < 0, (ii) 
metals with temperature-independent Pauli paramag- 
netic susceptibility xp > ^nd temperature-independent 
Landau diamagnetic susceptibility xl < well below the 
Fermi temperature T <^ Tp ~ Ep/kB, and (iii) magnetic 
materials with local moments with strongly temperature- 
dependent Curie paramagnetic susceptibility oc 1/T 
or Curie- Weiss law xiT) oc 1/{T — Tq) corresponding to 
ferromagnetic {Tq > 0) or antiferromagnetic {Tq < 0) 
interactions between the magnetic moments. When an 
electron configuration with J ~ (J: total angular mo- 
mentum) is realized for a magnetic ion, e.g., in the case 
of four d-electrons, Larmor diamagnetism and van-Vleck 
paramagnetism compete to determine the sign of x- 

Usually, the orbital motion of electrons leads to a dia- 
magnetic Xi BlS in the case of Larmor and Landau dia- 
magnetism. This is due to the cyclotron motion induced 
by the Lorentz force, which tends to reduce the action 
of the external magnetic field. When the band gap is 
reduced or fully closed, this diamagnetic susceptibility is 
considerably enhanced as in the cases of bismuth ||^ and 
graphene |^ . Especially in graphene, due to the presence 
of gapless two-dimensional Dirac fermions, x diverges as 
oc —1/T as T — > Despite this general tendency. 



there is no proof that the orbital motion always favors 
diamagnetism, in particular in the presence of strong 
spin-orbit interaction (SOI). Hence one cannot exclude 
the possibility that orbital paramagnetism may be real- 
ized in band insulators and metals. So far, however, this 
possibility has been considered by only a few theoreti- 
cal studies For example, Vignale Q showed that 
a two-dimensional (2D) electron gas in a periodic po- 
tential exhibits orbital paramagnetism when the Fermi 
energy is close to a saddle point of the band structure. 
Boiko and Rashba ||^ predicted orbital magnetism in the 
Rashba model, for which we will provide the first mate- 
rial realization and its further theoretical analysis below. 
For band insulators, it has so far been an open question 
whether orbital paramagnetism can be realized or not. 

In this paper, we address the mechanisms of orbital 
dia- and paramagnetism in semiconductors by studying 
BiTel, a narrow-gap polar semiconductor with a non- 
centrosymmetic layered structure. Due to the polarity 
of the bulk material and the presence of Bi atoms with 
large atomic SOI, BiTel exhibits an extraordinarily large 
bulk Rashba spin splitting (RSS), astonishingly reach- 
ing ~ 400 mcV at the bottom of the conduction band 
as observed in the recent angle-resolved photoemission 
spectroscopy experiment by Ishizaka et al. Further 
studies based on optical spectroscopy |]l0| , first-principles 
calculations and group theoretical analysis [|ll| have con- 
firmed that this giant RSS is indeed a bulk property of 
BiTel, and revealed that in addition to the bottom con- 
duction bands (BCB's), the top valence bands (TVB's) 
of BiTel are subject to a comparable RSS. These unique 
features make BiTel an ideal medium to explore different 
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aspects of orbital magnetism in a band insulator system. 
Below we theoretically predict and expericntally confirm 
that in BiTel both an enhanced orbital diamagnetism 
and an unconventional orbital paramagnetism can be re- 
alized by tuning the Fermi energy Ep. This intriguing 
behavior is shown to be a direct consequence of the RSS. 

As described in Refs. ^, |ll|, BiTel has a minimum 
band gap not at the Brillouin zone (BZ) center, but 
around the hexagonal face center of the BZ, referred to 
as point A, as shown in Figure |^(a) . Due to the strong 
covalency and ionicity of Bi-Te and Bi-I bonds, respec- 
tively, the BCB's are dominated by Bi-6p states, whereas 
the TVB's arc predominantly of Tc-5p character with a 
partial contribution from l-5p states. The huge SOI of Bi 
accompanied by a strong negative crystal field splitting 
of TVB's leads BCB's and TVB's to be symmetrically 
of the same character, and hence to be strongly coupled 
with each other via a quasi- two-dimensional RSS Hamil- 
tonian. Near point A, this can be described as 

HR = ^ + Xe,-{sx p), (1) 

where m* is the effective mass of the carriers, A is the 
Rashba parameter, is the unit vector in z direction, 
which is the direction of the potential gradient breaking 
the inversion symmetry, and s and p are the spin and 
momentum operators, respectively. Note in particular 
the crossing of the BCB's shown in Figure |l|(a), which 
is due to the Kramer's degeneracy at the time-reversal 
symmetry point A. As fitting parameters, m* and A can 
be tuned Q such that Hr can properly reproduce the 
electronic dispersion of BCB's around point A, especially 
the momentum offset of the conduction band minimum 
(CBM), IpI = Hko with ko = 0.052 A-\ and the Rashba 
energy Efi = 113 meV, where En is the energy of the 
conduction bands at their crossing point with respect to 
CBM. 

To calculate the orbital magnetic susceptibility x we 
employ the Fukuyama formula |l^ , 

x(T) = lio^^-^ksT'^'^TrlGiVr.GiVyGeVr.GiVy], 

e k 

(2) 

where Nj^ is the Avogadro constant, the number of 
carriers, Ge ~ [iu!£ + Ep — H]^^ the one-particle ther- 
mal Green's function with Matsubara frequencies uji = 
{2£ + l)7rfcBT, £ e Z, and Vi = dH/dh {i = x,y,z) 
the velocity operators. For the above 2D Rashba model. 
Figure |l|(b) shows the calculated dependence of x on the 
Fermi energy Ep. Near the band crossing point, where 
we take the origin of Ep, an enhanced diamagnetic sus- 
ceptibility is observed, which diverges as x c< —1/T as 
in the case of graphene For Ep » Er, x approaches 
— ^nioiMoe^/(127rm*), which is equivalent to the value of 
the Landau diamagnetic susceptibility for free electrons 
with mass m* in 2D (here ^moi is the molar surface area 
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FIG. 1: (Color online) (a) Band dispersion of BiTel near 
point A towards the H and L directions. The crossing point 
of the conduction bands is set at zero energy, (b) Orbital 
magnetic susceptibility obtained from a 2D Rashba model, 
which reproduces the dispersion of the BCB's in BiTel around 
point A. The main panel shows the Fermi energy dependence 
of X at r = 30 K. The arrow marks the value of the Landau 
diamagnetic susceptibility for free electrons with mass m* in 
2D (see the related discussion). The inset shows the temper- 
ature dependence of the diamagnetic peak value at Ep = 0, 
diverging as —1/T. 



of one BiTel layer). In addition, we find the paramag- 
netic orbital susceptibility when Ep is below the band 
crossing point, which is a unique feature attributed to 
the Rashba spin splitting . 

In the zero temperature limit, by evaluating the fre- 
quency sum in the Fukuyama formula we have obtained 
an analytic formula for the orbital susceptibility of the 
2D Rashba model, 

e2 1 f 8fc2-3fc^ 

X(T = 0) = ,0^..-^ j ^^dfc. 

E-{k}<EF<E+{k) 

Here k = |k| denotes the radial momentum variable, 
while the angular integration has already been per- 
formed; E±{k) = h^k^/2'm* ± Afc are the two branches 
of the Rashba dispersion and /co ^ m*X/lT? the mini- 
mum of the lower branch. For small band fillings, i.e., 
Ep = —En, the momentum integral is limited to the re- 
gion where k = ko, hence x is positive and approximately 
given by 

X(T = 0) - /io A,noi ^ tJ- v/1 + Ep/Er. (4) 

On the other hand, if the Fermi energy is near the band 
crossing, i.e., Ep = 0, the region of small k gives the 
dominant contribution, which gives rise to the negative 
divergence of the susceptibility as x oc —l/\Ep\. 

As a supporting argument for the orbital paramag- 
netism at low band fillings, we have shown in the sup- 
plemental material that the energy of the lowest Landau 
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level of the 2D Rashba model in a perpendicular magnetic 
field B is always lower than the minimum energy —En in 
the absence of B. Accordingly, for sufficiently low filling 
of the bands, by applying the magnetic field the total en- 
ergy of the system is lowered, a sign for the occurrence of 
paramagnetism. Since the Landau-Peierls formula [ p^ , 
which takes into account only intraband contributions, 
gives always a negative value for the orbital magnetic 
susceptibility, the positive x can be attributed to the in- 
terband contributions, which are allowed due to the RSS 
of BCB's. Also note that this behavior is different from 
that in the 2D Dirac model describing graphene, which 
shows only a diamagnetic orbital susceptibility. 

To further corroborate these results as well as to ex- 
plore the effects of all energy bands, we have performed 
a detailed analysis based on first-principles density func- 
tional theory (DFT). The orbital magnetic susceptibility 
of BiTel has been computed using the Fukuyama for- 
mula (^, with the velocity operator Vi {i = x,y,z) rep- 
resented in matrix form as 



'Ji,,un = {n\vi{k)\m) = ^(ii 



dHjk) 
dki 



(5) 



where |n) corresponds to the nth eigenstate of i?(k). The 
numerical computation of Ui.nm has been done using the 
approach of Ref. whereby i?(k) is described in the 
basis of so-called maximally localized Wannier functions 
(MLWF's) 0. To fully take into account the contribu- 
tions of all possible energy states, a set of 18 MLWF's, 
spanning the 12 highest valence bands and 6 lowest con- 
duction bands, has been generated with the wannierQO 
code 1 18 by postprocessing the relativistic DFT cal- 
culations, initially done using Perdew-Burke-Ernzerhof 
exchange-correlation functional ||2^ and the augmented 
plane wave plus local orbitals method as implemented 
in the wien2k package |2^. For these calculations, the 
muffin tin radii are set to Rmt = 2.5 Bohr for all the 
atoms and the maximum modulus for the reciprocal vec- 
tors Kmax is chosen such that Rmt K max = 7.0. The BZ 
is sampled using a 20 x 20 x 20 k-mesh and a full structural 
optimization is performed until the magnitude of force on 
each ion is less than 0.1 mRy/Bohr (see Ref. 11 for more 
details). For the evaluation of the momentum integral in 
Eq. (H) we have used a 500 x 500 x 50 k-mesh to sam- 
ple BZ. We have also calculated the Pauli paramagnetic 
susceptibility xpi leaving the orbital degree of freedom 
unaffected. Namely we have added only the Zeeman cou- 
pling —^Bcr ■ B to the Hamiltonian and calculated the 
spin magnetization to obtain xp- However, the contribu- 
tion of is much smaller than the orbital x s-s described 
below. 

Figure ^ shows the Ep dependence of x + Xp obtained 
by the first-principles calculations for BiTel taking into 
account the dispersion along kz direction, where the ori- 
gin of Ep is taken at the crossing point of BCB's at 
point A. The Pauli contribution is an order of magni- 



tude smaller than the orbital x, at least in the considered 
range of -0.15 eY < Ep < 0.15 eV. This is due to the 
small effective mass of BCB carriers within the hexagonal 
face of the BZ (hereafter referred to as ab plane), m* ~ 

0. 2m |l^. This is analogous to the case of bismuth |^ 
and can be understood by considering a spin-degenerate 
two-dimensional parabolic band E = h^k^ /{2m*). Then 
the orbital and Pauli susceptibility are given by x = 
-{l/i)^ioNAfJ'*B'^/Ep and xp = ^J,oNA^J'B'^ / Ep, respec- 
tively, where fig ~ eh/ {2m*), fiB ~ eh/ {2m), and 
Ep = h'^kp/{2m*) with the Fermi wave vector kp. Since 
TO* appears in Ep, but not in the spin-interaction term, 

1. e. in /xs, it turns out that xp = — 3(rn*/TO)^x — ~0.1x. 
Therefore we can neglect xp focus on x below. 

For B parallel to e^, near the band crossing point a 
large temperature-dependent diamagnetic susceptibility 
is observed, while the susceptibility turns to be param- 
agnetic and temperature-independent when Ep is away 
from it. In contrast, if B is applied perpendicular to the 
z direction, x is found to be always positive and almost 
insensitive to Ep. The calculated results are in good 
agreement with the experiment, as shown in Figure ^ 
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FIG. 2: (Color online) Magnetic susceptibility of BiTel as a 
function of Fermi energy Ep, at two temperatures T = 300 K 
(upper panel) and T = 30 K (lower panel). Solid lines repre- 
sent the calculated susceptibility including the Pauli contri- 
bution, and points show the corresponding experimental data 
after subtracting the core contributions. Ep = corresponds 
to the crossing point of the Rashba-split conduction bands at 
point A. 
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For the experimental measurement of x, single crys- 
tals of BiTel with various charge carrier densities were 
grown by a Bridgman method or chemical vapor trans- 
port. The obtained crystals have cleaved (001) planes. 
The charge carrier density n of each sample was deter- 
mined by the Hall resistivity measurement, and the cor- 
responding Fermi energy Ep was deduced from the calcu- 
lated relationship between n and Ep in the 18-band tight- 
binding model, as shown in Figure 3(a) of Ref. |lO[ The 
magnetization measurement was performed with use of 
the vibrating-samplc magnetometer. The sample holder 
for this was made of a plastic straw tube, whose tiny 
diamagnetic signal was accurately subtracted. 

To extract the orbital magnetic susceptibility from the 
experimental data, one needs to subtract the Larmor sus- 
ceptibly XL originating from the individual ionic cores in 
BiTel. Since in this material the ionic states of Bi, Te 
and I ions are 3-t-, 2- and 1-, respectively, one can assume 
that they all have closed shell ionic configurations. Ac- 
cordingly, for each ion xl can be easily estimated as [|] , 
XL = -Aio^A(e^/6TO) X;i=„,ij 2'j(r2)j, where is the 
number of electrons occupying state i = {n,l,j}. Per- 
forming a set of all-electron DFT calculations using the 
LDAl.x program [Q, the spread functions (r^)i have 
been individually calculated for all the occupied states 
of Bi'^+, Te^~ and (see Table 1 in the supplemen- 
tal material). Based on these calculations, the respective 
XL of Bi3+, Te^^ and I^" are found to be -32.5, -23.3 
and -21.1 (all in units of 10~^ emu/mol). We have ac- 
cordingly subtracted these values from our experimental 
data to deduce the orbital magnetic susceptibility. 

In Figure |^, we show the temperature dependence 
of the orbital magnetic susceptibility for B || at 
two characteristic Fermi energies, one near the band 
crossing point {Ep = 0) and the other away from it 
{Ep ~ 70 meV). In the first case, a strong tempera- 
ture dependence and sign change of x from paramag- 
netism to diamagnetism is observed as the temperature 
is lowered, while in the latter case the susceptibility re- 
mains paramagnetic and almost independent of tempera- 
ture. This behavior is clearly understood in terms of the 
temperature-dependent enhanced diamagnetic contribu- 
tion near Ep = 0, which adds to the almost temperature- 
independent paramagnetic interband contribution. All 
these features are in excellent agreement with the exper- 
imental curves measured at two different samples with 
corresponding carrier densities, as shown in Fig ^ 

Up to now, we have focused on the experimentally ac- 
cessible region — 0.15 eV < i?F < 0.15 eV. It is of interest 
to study X in extended region to reveal the mechanism 
of the enhanced orbital magnetism. Figure ^a) shows 
the calculated orbital magnetic susceptibility as a func- 
tion of Ep over the whole energy range of the 18-band 
model at T 300 K. Large positive and negative values 
are observed, which arc even further enhanced if we con- 
sider only the contribution from the ab plane of the BZ 
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FIG. 3: (Color online) Temperature dependence of the orbital 
magnetic susceptibility of BiTel for B || e^, obtained theo- 
retically (left panel) and experimentally (right panel) at two 
characteristic carrier densities. For n — 3.4 x lO^^cm"'^ the 
Fermi energy Ep is close to the band crossing point, whereas 
for n = 6.0 x lO^^cm"^ it is ~ 70 meV above it. 



as in Figure ^b). It turns out that these large dia- and 
paramagnetic peaks are caused by (anti)crossings in the 
band structure of BiTel, which can be described as 2D 
tilted Dirac cones. A detailed discussion of the appear- 
ance of such tilted Dirac cones in the band structure of 
BiTel and their correspondence to the observed dia- and 
paramagnetic peaks in Figure ^ is given in the supple- 
mental material. There we also provide a simple model 
which generically describes the dispersion of tilted band 
crossings and explains the diverging paramagnetic orbital 
susceptibility as the temperature T goes to zero. 

In summary, we have studied the orbital magnetism fo- 
cusing on its enhancement due to interband effects. As a 
concrete example we have considered BiTel with its gi- 
ant Rashba splitting to see how the spin-orbit interac- 
tion affects the orbital magnetism. In addition to the 
temperature-dependent large diamagnetic susceptibility 
near the band crossing point analogous to graphene, we 
have found two mechanisms for the enhanced orbital 
paramagnetism. One is the Rashba splitting with the 
Fermi energy below the band crossing, and the other is 
the tilted band crossing. 
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FIG. 4: (Color online) Theoretical results for the orbital mag- 
netic susceptibility in the entire energy range of the 18-band 
model of BiTel at T = 300 K. For (a) we have performed 
in Eq. (^) the k integral over the whole three-dimensional 
BriouUin zone, whereas for (b) we have calculated only the 
contribution from the ab plane, i.e., we have assumed the 
integrand to be indepent of and taken only the value at 

fcz = Tt/c. 



